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ABSTRACT 

Number counts of galaxies are re-analyzed using a semi-analytic model (SAM) of galaxy formation based on 
the hierarchical clustering scenario. We have determined the astrophysical parameters in the SAM that reproduce 
observations of nearby galaxies, and used them to predict the number counts and redshifts of faint galaxies for 
three cosmological models for (1) the standard cold dark matter (CDM) universe, (2) a low-density flat universe 
with nonzero cosmological constant, and (3) a low-density open universe with zero cosmological constant. The 
novelty of our SAM analysis is the inclusion of selection effects arising from the cosmological dimming of surface 
brightness of high-redshift galaxies, and also from the absorption of visible light by internal dust and intergalactic 
H I clouds. Contrary to previous SAM analyses which do not take into account such selection effects, we find, 
from comparison with observed counts and redshifts of faint galaxies in the Hubble Deep Field (HDF), that the 
standard CDM universe is not preferred, and a low-density universe either with or without cosmological constant 
is favorable, as suggested by other recent studies. Moreover, we find that a simple prescription for the time scale 
of star formation (SF), being proportional to the dynamical time scale of the formation of the galactic disk, is 
unable to reproduce the observed number- redshift relation for HDF galaxies, and that the SF time scale should be 
nearly independent of redshift, as suggested by other SAM analyses for the formation of quasars and the evolution 
of damped Ly-a systems. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: formation - large-scale structure of 
universe 



1. INTRODUCTION 

In the field of observational cosmology, counts of the number 
of faint galaxies in a given area of sky provide one of the most 
fundamental observables from which the cosmological parame- 
ters, and hence the geometry of the Universe, can be determined 
(e.g., Peebles 1993). The predicted number counts of galaxies 
as a function of apparent magnitude and redshift, which should 
be compared with the data, is obtained from summing up the 
derived luminosity functions over all redshifts and morphologi- 
cal types, then multiplied by a redshift-dependent cosmological 
volume element. While the local luminosity function of galax- 
ies of individual types is known from redshift surveys, the lu- 
minosity function at any high redshift still has to be deduced 
with the help of morphological-type-dependent evolution mod- 
els of galaxies. Accordingly, the predicted number counts rely 
directly on how well the evolution of galaxies is modeled from 
their formation until the present (e.g., Yoshii & Takahara 1988). 

With the usual assumption of monolithic collapse, the wind 
model for elliptical galaxies and the infall model for spiral 
galaxies are able to reproduce many of the observed proper- 
ties of nearby galaxies, and provide a strong theoretical tool for 
understanding their evolution (Arimoto & Yoshii 1986, 1987; 
Arimoto, Yoshii, & Takahara 1991). However, when using these 
traditional evolution models, it has been claimed that the stan- 
dard cold dark matter (CDM) universe, or the Einstein-de Sit- 
ter (EdS) universe, is not reconciled with the observed high 
counts of faint galaxies, and that a low-density universe is pre- 
ferred (Yoshii & Takahara 1988; Yoshii & Peterson 1991, 1995; 



Yoshii 1993). Most recently, Totani & Yoshii (2000, hereafter 
TY00) compared their predictions against the number counts 
observed to the faint limits in the Hubble Deep Field (HDF, 
Williams et al. 1996) by taking into account various selection 
effects for the first time. Allowing for the possibility of num- 
ber evolution of galaxies in a phenomenological way, and af- 
ter a comprehensive check of systematic model uncertainties, 
they strengthened previous claims and further demonstrated 
that mild or negligible merging of high-redshift galaxies in a 
low-density flat universe with nonzero cosmological constant 
is a likely solution to simultaneously reproduce the observed 
high counts and redshift distribution of faint galaxies. On the 
other hand, the EdS universe is in serious contradiction with the 
data even if a strong number evolution is invoked. 

However, there is a growing evidence from recent observa- 
tions of the large-scale structure of the universe that gravitation- 
ally bound objects, such as galaxy clusters, are formed through 
continuous mergers of dark halos. Some authors have tried to 
construct a galaxy formation model based on this scenario of hi- 
erarchical clustering in a CDM universe, referred to as a semi- 
analytic model (SAM) of galaxy formation (e.g., Kauffmann, 
White, & Guiderdoni 1993; Cole et al. 1994). In fact, SAMs 
successfully reproduced a variety of observed features of local 
galaxies such as their luminosity function, color distribution, 
and so on. 

The number counts of faint galaxies have also been analyzed 
using SAMs (Cole et al. 1994; Kauffmann, Guiderdoni & 
White 1994; Heyl et al. 1995; Baugh, Cole & Frenk 1996). 
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The Durham group claimed that their model agrees with the 
observed number counts in the EdS universe rather than in a 
low-density universe. The Munich group reached a similar con- 
clusion by halting the process of star formation in smaller halos 
with circular velocities less than 100 km s" 1 . These conclusions 
by both groups are apparently at variance with those obtained 
using the traditional evolution models of galaxies. This appar- 
ent discrepancy needs to be immediately resolved. 

In this paper we examine the importance of selection effects 
in observations of faint galaxies which have been ignored in 
previous SAM analyses of galaxy number counts, such as the 
cosmological dimming of surface brightness and the absorption 
of emitted light by internal dust and intergalactic H I clouds. 
Taking into account these effects in the same way as in TY00, 
our SAM analysis obtains, for the first time, the predicted num- 
ber counts that can be consistently compared with the HDF 
counts to faint magnitude limits. 

This paper is outlined as follows. In §2 we briefly describe 
our SAM. In §3 we constrain the astrophysical parameters in 
our SAM analysis against local observations. In §4 we compare 
theoretical number counts and redshift distributions of faint 
galaxies with the data, and discuss the range of uncertainties 
in our calculations of galaxy counts. In §5 we provide a sum- 
mary and discussion. 

2. MODEL 

The SAM we employ involves known physical processes 
connected with the process of galaxy formation. It is therefore 
straightforward to understand how galaxies form and evolve 
within the context of this model. In the CDM universe, dark 
matter halos cluster gravitationally and are merged in a manner 
that depends on the adopted power spectrum of initial density 
fluctuations. In each of the merged dark halos, radiative gas 
cooling, star formation, and gas reheating by supernovae oc- 
cur. The cooled dense gas and stars constitute galaxies. These 
galaxies sometimes merge together in a common dark halo and 
more massive galaxies form. 

For the purpose of comparison with observation, we use a 
stellar population synthesis approach, from which the luminosi- 
ties and colors of model galaxies are calculated. The SAM well 
explains the local luminosity function of galaxies, the color dis- 
tribution, and so on. Our present SAM analysis obtains essen- 
tially the same results of previous SAM analyses, with minor 
differences in a number of details. In this section we only 
briefly describe our model; its full description will be given in 
Nagashima & Gouda (2001, in preparation). 

2. 1 . Scheme of Galaxy Formation 

First, based on the method of Somerville & Kolatt (1999), 
we determine the merging history of dark matter halos by 
extending the Press-Schechter formalism (Press & Schechter 
1974; Bower 1991; Bond et al. 1991; Lacey & Cole 1993). 
We adopt the power spectrum for a specific cosmology from 
Bardeen et al. (1986), and assume a halo with circular veloc- 
ity Vcirc <40km s" 1 as a diffuse accretion entity. The evolution 
of the baryonic component is followed until the output redshift 
with the redshift interval of Az = 0.06(1 +z), corresponding to 
the dynamical time scale of halos which collapse at that time. 
In order to minimize the artificial effect of dividing the his- 
tory of galaxy formation into discrete redshift intervals, we fix 
Az = 0.06 until just prior to the output redshift. This manipula- 
tion is important especially at high redshift. In our analysis, the 



highest output redshift is z — 9 with the interval of the output 
redshift ~0.05-0.4. 

If a dark matter halo has no progenitor halos, the mass frac- 
tion of the gas is given by ilb/^o > where fio is the den- 
sity parameter for the total baryonic and non-baryonic com- 
ponents, and fib = 0.015/T 2 is the baryonic density parameter 
constrained by primordial nucleosynthesis calculations (e.g., 
Suzuki, Yoshii, & Beers 2000). In this work h represents the 
Hubble parameter given by h = Ho/ 100 km s _1 Mpc _1 . 

When a dark matter halo collapses, the gas in the halo is 
shock-heated to the virial temperature of the halo. We refer to 
this heated gas as the hot gas. At the same time, the gas in dense 
regions of the halo is cooled due to efficient radiative cooling. 
We call this cooled gas the cold gas. Assuming an isothermal 
density distribution of the entire halo and using the metallicity- 
dependent cooling function by Sutherland & Dopita (1991), we 
calculate the amount of cold gas which eventually falls onto a 
central galaxy in the halo. In order to avoid the formation of un- 
physically large galaxies, the above cooling process is applied 
only to halos with V c i rc <500 km s" 1 in the standard CDM and 
400 km s" 1 in low-density universes. The physical reason for 
this restriction is not clear, but Cole et al. (2000a) recently re- 
ported that the formation of enormous galaxies is hindered in 
the isothermal halo having a large core. Nevertheless, in this 
paper we adopt a simple isothermal distribution and prevent the 
formation of so-called "monster galaxies" by hand. 

Stars are formed from the cold gas at a rate of M. t = M co id / t* , 
where M co id is the mass of cold gas and r* is the time scale of 
star formation. We assume that r* is independent of z, but de- 
pendent on y c i rc as follows: 
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This form is referred to as the "Durham model" by 
Somerville & Primack (1999). It should be noted that Cole et 
al. (2000a) modified their original form (equation 1) as 

Vcirc Y' ( T ^( Z ) \ (2 ) 

,200km s" 1 / Udyn(0)y' 
by multiplying a factor proportional to the dynamical time scale 
in the galactic disk. This form is similar to "Munich model" 
(Kauffman, White & Guiderdoni 1993). 

For convenience, the cases in equations (1) and (2) are here- 
after referred to as "constant star formation (CSF)" and "dy- 
namical star formation (DSF)," respectively. The free param- 
eters of t° and a* in CSF and DSF are fixed by matching the 
observed mass fraction of cold gas in neutral form in the disks 
of spiral galaxies. Cole et al. (2000a) stressed that the V C i rc - 
dependence is needed to reproduce the observed ratio of gas 
mass relative to the B-band luminosity of a galaxy. The effect 
of introducing the Td yn -dependence in DSF will be discussed in 
§4.2.1. 

In our SAM, stars with masses larger than 10M Q explode 
as Type II supernovae (SNe) and heat up the surrounding cold 
gas. This SN feedback reheats the cold gas to hot gas at rate of 
Reheat = Areola /^reheat, where the time scale of reheating is given 
by 



^reheat : 



^circ 



(3) 



The free parameters of Vhot and «hot are determined by 
matching the local luminosity function of galaxies. 
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With these M* and M re heat thus determined, we obtain the 
masses of hot gas, cold gas, and disk stars as a function of time 
during the evolution of galaxies. Chemical enrichment is also 
taken into account adopting heavy-element yield ofy = 0.038 = 
2Z©, but changing this value of y has a minimal effect on the 
results described below. 

When two or more progenitor halos have merged, the newly 
formed larger halo should contain at least two or more galaxies 
which had originally resided in the individual progenitor halos. 
By definition, we identify the central galaxy in the new com- 
mon halo with the central galaxy contained in the most massive 
of the progenitor halos. Other galaxies are regarded as satellite 
galaxies. 

These satellites merge by either dynamical friction or ran- 
dom collision. The time scale of merging by dynamical friction 
is given by 

260 / R \ 2 f V arc \ ( M sat y 1 ^ 

where R and y c ; rc are the radius and the circular velocity of 
the new common halo, respectively, In A c is the Coulomb loga- 
rithm, and M sat is the mass of the satellite galaxies including the 
mass of dark matter (Binney & Tremain 1987). When the time 
elapsed after a galaxy becomes a satellite exceeds tm c . a satel- 
lite galaxy infalls onto the central galaxy. On the other hand, 
the mean free time scale of random collision is given by 

_ 500 / R \ 3 f r gd V 2 
Tco11 N 2 \MpcJ ^0.12Mp C y 

x (_£LY7_f!-^yGyr > (5) 
V 100km s" 1 / V 300km s" 1 / 

where N is the number of satellite galaxies, r ga i is their ra- 
dius, and (Thdo and cr ga i are the ID velocity dispersions of the 
common halo and satellite galaxies, respectively (Makino & 
Hut 1997). For simplicity, the satellite radius r ga i is set to be 
one tenth of the virial radius of a progenitor halo to which the 
satellite once belonged as a central galaxy. With a probability 
Af/r co u, where Af is the time step corresponding to the redshift 
interval Az, a satellite galaxy merges another satellite picked 
out randomly. This process was first introduced in a SAM by 
Somerville & Primack (1999). 

Consider the case that two galaxies of masses m\ and m^iy- 
mi ) merge together. If the mass ratio f = m\ /ni2 is larger than a 
certain critical value of /bulge. we assume that a starburst occurs 
and all the cold gas turns into hot gas, which fills in the halo, 
and the stars populate the bulge of a new galaxy. On the other 
hand, if / < /bulge, no starburst occurs and a smaller galaxy is 
simply absorbed into the disk of a larger galaxy. These pro- 
cesses are repeated until the output redshift. 

Given the SF rate as a function of time or redshift, the ab- 
solute luminosity and colors of individual galaxies are calcu- 
lated using a population synthesis code by Kodama & Arimoto 
(1997). The stellar metallicity grids in the code cover a range 
of Z* =0.0001-0.05. The initial stellar mass function (IMF) that 
we adopt is the power-law IMF of Salpeter form with lower 
and upper mass limits of 0. 1M Q and 60M Q , respectively. Since 
our knowledge of the lower mass limit is incomplete, there is 
the possibility that many brown dwarf-like objects are formed. 
Therefore, following Cole et al. (1994), we introduce a param- 
eter defined as T = (Mi um +M B D)/Mi um , where M\ um is the total 
mass of luminous stars with m > 0.1M Q and M BD is that of 
invisible brown dwarfs. 



To account for extinction by internal dust we adopt a simple 
model by Wang & Heckman (1996) in which the optical depth 
in B-band is related to the luminosity as t b = 0.8(L B /1.3 x 
10 10 L Q ) - 5 . Optical depths in other bands are calculated by us- 
ing the Galactic extinction curve, and the dust distribution in 
disks is assumed to be the slab model considered by Somerville 
& Primack (1999). It is not trivial that such an empirical dust 
model can be extraporated to very bright galaxies. However, 
since the extinction is typically Ab ~ 1 mag at the bright-end 
of the local luminosity function, our results especially in the 
/-band are not significantly affected by the details of the dust 
model. 

Emitted light from distant galaxies is absorbed by Lyman 
lines and Lyman continuum in intervening intergalactic H I 
clouds. The redshift at which this effect becomes important is 
different for different photometric passbands, because the ab- 
sorption occurs around 1000 A in the rest frame of the clouds. 
Fig. Al shows the expected extinction as a function of red- 
shift for four passband filters of the HDF (Yoshii & Peterson 
1994). This effect has been included into a SAM by Baugh et 
al. (1998) in order to pick out the Lyman-break galaxies at high 
redshift by color selection criteria. 

We classify galaxies into different morphological types ac- 
cording to the B-band bulge-to-disk luminosity ratio B/D. In 
this paper, following Simien & de Vaucouleurs (1986), galax- 
ies with B/D > 1.52, 0.68 <B/D< 1.52, and B/D < 0.68 are 
classified as ellipticals, SOs, and spirals, respectively. Kauff- 
mann et al. (1993) and Baugh, Cole & Frenk (1996) showed 
that this method of type classification well reproduces the ob- 
served type mix. 

The above procedure is a standard one in the SAM of galaxy 
formation. In the next subsection, in order to investigate the 
properties of high-z galaxies properly, we introduce two impor- 
tant effects into our SAM analysis. 

2.2. New ingredients of the model: selection effects 

We judge whether the surface brightness of the galaxies in 
our SAM is above the detection threshold of the HDF obser- 
vations. The intrinsic size of spiral or disk-dominated galax- 
ies is estimated by adopting the dimensionless spin parameter 
Ah = 0.05, and by conserving the specific angular momentum 
during the gas cooling. On the other hand, the intrinsic size 
of early-type galaxies is estimated from their virial radii, ad- 
justed by a scaling parameter fi, to match with the observed 
size of early-type galaxies. Note that the selection effects on 
the local luminosity function have been considered by Cole et 
al. (2000a) in a simple way in which the isophotal magnitude 
within 25 mag arcsec" 2 was used. 

The selection effects in predicting the HST number counts in 
our SAM analysis are evaluated as follows. Using the intrinsic 
size of model galaxies as obtained above, and adopting a Gaus- 
sian point-spread function f(x) for the HST observations, the 
surface brightness profile g(x) in the observer frame is given by 

g(\x\)=J^'f(\x'-x\)g(\x'\), (6) 

where |x| = x = r/r e is the normalized radius away from the 
galaxy center relative to the effective radius r e , and g(x) = 
exp(— a n x l ' n ) is the intrinsic surface brightness profile. We 
adopt n = 1 for spirals and n = 4 for ellipticals. The coeffi- 
cient a n is given by a\ = 1.68 and = 7.67. Then, for each 
of our model galaxies, we determine their surface brightness 
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profile S(8), where 8 = xr e /dA and d\ is the angular-diameter 
distance. We note that model galaxies with surface brightness 
brighter than the threshold 5 t h and with an isophotal diameter 
larger than the minimum diameter D min are actually detected 
as galaxies (Yoshii 1993). In order to be consistent with the 
HDF observations, we use the isophotal magnitude scheme, 
5th = 27.5 mag arcsec" 2 in Veoe, 5th = 27.0 mag arcsec" 2 in 
h\4, U300 and B450, and D min ~0.2 arcsec. More details are de- 
scribed in TY00. 

3. THE SETTING OF PARAMETERS IN OUR SAM ANALYSIS 

We consider the predicted number counts in four models - 
SC, OC, LC, and LD (Table Al). The first three models are for 
CSF (equation 1), and the fourth model is for DSF (equation 2). 
The capitals S, O, and L refer to the standard CDM universe, a 
low-density open universe and a low-density flat universe with 
non-zero cosmological constant (A), respectively. 

The cosmological parameters (fio, ^A, h, a$) are tabulated 
in Table Al. For all the models the baryon density parameter 
fib = 0.015/z" 2 is used in common. For the low-density open 
and flat universes the value of cr% is determined from observed 
cluster abundances (Eke, Cole & Frenk 1996). 

The astrophysical parameters (Vh ot , Ofhot, t®, a*, f b , /bulge, T) 
are constrained from local observations. However, since the pa- 
rameter of /bulge, among others, do not affect our result, we set 
/bulge = 0.2 for the standard CDM universe, and /bulge = 0.5 for 
the low-density open and flat universes. Other parameters are 
discussed below. 

3.1. The Local Luminosity Function of Galaxies 

The SN feedback-related parameters of Vh t and ahot deter- 
mine the location of the knee of the luminosity function and the 
faint-end slope, respectively. It should be noted that the mass 
fraction T of invisible stars determines the magnitude scale of 
galaxies, so that changing T moves the luminosity function hor- 
izontally without changing its overall shape. Therefore, cou- 
pled with Vhot, Y determines the bright portion of the luminosity 
function. 

Fig. A2 shows theoretical results represented by thick lines 
for the four models tabulated in Table Al. Symbols with er- 
rorbars indicate observational results from the B-band redshift 
surveys such as APM (Loveday et al. 1992), ESP (Zucca et al. 
1997), Durham/UKST (Ratcliffe et al. 1997), 2dF (Folkes et al. 
1999) and SDSS (Blanton et al. 2000), and from the #-band 
redshift surveys (Mobasher et al. 1993; Gardner et al. 1997; 
2MASS, Cole et al. 2000b). It is evident that while most of B- 
band redshift surveys give a rather flat slope in the faint end, the 
ESP and SDSS surveys give an extreme case showing a much 
steeper slope. Note that the SDSS luminosity function shown 
in Fig. A2 is that with the same detection limit as employed in 
the 2dF survey. 

In this paper we have chosen the values of SN feedback- 
related parameters so as to reproduce a flat faint-end slope, but 
in §4.2.3 we will investigate the effect of using the steepest ESP 
slope in predicting the number counts of faint galaxies. 

3.2. The Mass Fraction of Cold Gas in Spiral Galaxies 

The SF rate-related parameters of r» and a* determine the 
overall mass fraction of cold gas in galaxies with given lumi- 
nosity and its luminosity dependence, respectively. In some 
previous SAM analyses, the mass fraction of gas in the Milky 
Way is exclusively used to fix r° for all other spiral galaxies. 



Then, when a* = 0, Cole et al. (2000a) found that dwarf-size 
galaxies have too little gas to be consistent with observations. 
According to Cole et al. (2000a), we here constrain both r° 
and a* in a combined manner to reproduce the mass fraction of 
cold gas in galaxies spanning a wide range of luminosity. 

Fig. A3 shows the ratio of cold gas mass relative to B-band 
luminosity of spiral galaxies as a function of their luminosity. 
Solid curves show theoretical results for the four models tabu- 
lated in Table Al. We here assume that 75% of the cold gas in 
these models is comprised of hydrogen, i.e., Mm = 0.75M co id. 
Filled diamonds with errorbars indicate the H I data taken from 
Huchtmeier & Richter (1988). 

It should be noted that use of t® for either CSF (equation 
1) or DSF (equation 2) hardly affects the resulting luminosity 
function, as seen from the difference between LC and LD in 
Fig. A2. However, a* depends on the strength of the SN feed- 
back, and the strength should be different for different universe 
models (see Table Al). Therefore, in the case of a*, we can 
only fix its value after the SN feedback-related parameters are 
constrained in a specified universe model. 

3.3. The Intrinsic Sizes of Galaxies 

The scaling parameter /, is determined by matching the ef- 
fective radius of early-type galaxies, while the actual In- 
dependence in /, is ignored for simplicity. We discuss the effect 
of this assumption in §4.2.2. 

Knowledge of the intrinsic size of galaxies is a necessary 
quantity in our SAM analysis, because it is the surface bright- 
ness of galaxies, not their luminosity, that is relevant to evaluat- 
ing the selection effects for model galaxies. For spiral galaxies, 
the key process is the conservation of specific angular momen- 
tum during the cooling of hot gas. The effective disk radius r e is 
given by r e = (1.68/V2)AH?"i, where r\ is an initial radius of the 
progenitor gas sphere (Fall 1983). However, the definition of 
disk size during the cooling phase is not readily apparent from 
this formula. Furthermore, the angular momentum transfer dur- 
ing mergers of galaxies in a merged halo is very complicated. 
Without entering into all these complexities, in this paper we 
re-estimate the disk size by the above equation when the disk 
mass increases twice. This approach is simple but reproduces 
the observations rather well. 

Fig. A4 shows the effective disk radii of spiral galaxies as a 
function of their luminosity. Thick solid lines show the theoret- 
ical results for the four models tabulated in Table Al. Dotted 
lines indicate the best-fit relation to the observational data given 
by Impey et al. (1996). All models provide reasonable disk 
sizes for dwarf galaxies. Slight deviations from the observa- 
tions of bright galaxies may have occurred partly because dust 
obscuration becomes effective for such galaxies (cf. §2.2), and 
partly because our method of defining the disk size is too sim- 
ple. We discuss the effect of changing the disk size in §4.2.2. 

In this paper, the effective radius of elliptical galaxies is es- 
timated from scaling the virial radius by a parameter /,, i.e., 
r e = fbGM\,/V? rc , where Mb is the mass of stars and cold gas, 
and /, is fixed to reproduce the effective radius of an galaxy. 
In Fig. A4, thick dashed lines indicate the effective radius 
of elliptical galaxies for the four models tabulated in Table 
Al. Dashed and dot-dashed lines show the best-fit relations for 
dwarf and compact ellipticals, respectively, based on the obser- 
vational data given by Bender et al. (1992). Unfortunately, our 
SAM analysis cannot reproduce two distinct sequences simul- 
taneously; the theoretical r e - L B relation becomes coincident 
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with the bright portion of the giant-dwarf sequence, whereas 
the same relation changes its slope becoming coincident with 
the faint portion of giant-compact sequence. We find that this 
changeover magnitude is determined by the strength of SN 
feedback (see Appendix A), and that the extreme manipulation 
that reproduce the dwarf sequence affects the number counts of 
model galaxies only weakly. 

4. RESULTS 
4. 1 . Galaxy Number Counts 

Fig. A5 shows the number counts of galaxies as a function 
of apparent isophotal magnitude for the SC, OC, and LC mod- 
els in the HST UBVI bands in the AB magnitude system. The 
thick lines are the theoretical predictions, based on the HST 
observational conditions, including the selection effects from 
the cosmological dimming of surface brightness and also from 
the absorption of visible light by internal dust and intergalactic 
H I clouds. The thin lines are the predictions ignoring the se- 
lection effects except for the effect of dust absorption, because 
the dust absorption is already taken into account in reproducing 
the local luminosity function. Open circles with errorbars show 
the HDF data (Williams et al. 1996), and other symbols show 
ground-based data after transformation to AB magnitudes to be 
consistent with the HST. 

In order to discriminate favorable universe models, the pre- 
dictions shown by the thick lines should be compared directly 
with the observational data. It is evident from this figure that the 
SC model falls short of the HST data. We note that the discrep- 
ancy between the SC model and the data is seen at B450 > 25, 
even if we do not consider the selection effects. On the other 
hand, both of the LC and OC models agree well with the data, 
owing to the inclusion of selection effects. We will discuss the 
uncertainty of estimating these effects in the next section. 

Fig. A6 shows the redshift distribution of galaxies. The thick 
lines denote the SC, OC and LC models, including all the se- 
lection effects as in Fig. A5. The histogram indicates the ob- 
served redshift distribution based on photometric redshifts of 
HDF galaxies estimated by Furusawa et al. (2000), in which 
they improved the method of redshift estimation compared to 
the method by Fernandez-Soto, Lanzetta & Yahil (1999). 

In the LC model, both the peak height and the distribution 
towards higher redshift agree well with the observation. How- 
ever, the peak height of the SC model falls significantly short 
of that observed at z — 1 - 1 .5. This is a direct reflection of the 
lack of galaxies in the theoretical number-magnitude relation 
for the SC model (see Fig.A5). In the OC model, the relative 
number between the peak and high-z tail at z > 3 is smaller 
than that observed and is inconsistent with the data. This is a 
direct reflection of the z-dependence of the comoving volume 
element dV /dz, which is shown in the middle panel of Fig.A6. 
While in the LC model the matter density dominates over the 
cosmological constant at high redshift, which leads to a similar 
z-dependence of dV /dz to that in the EdS universe, the nega- 
tive curvature effect makes dV /dz decline more slowly at such 
high redshift in the OC model. Thus the redshift distribution in 
the OC model declines toward higher redshift more slowly than 
that in the LC model. In order for the OC model to agree with 
the data, we would need to halt the star formation at z > 3 by 
hand, for unknown reasons. 

We conclude that as far as that the astrophysical parameters 
in our SAM analysis are constrained by local luminosity func- 
tion of galaxies, the standard CDM universe does not agree with 



both the observed number-magnitude and number-redshift re- 
lations, and that the A-dominated flat universe is best able to 
reproduce these relations simultaneously. 

4.2. Uncertainties in the SAM counts 

In this section we discuss uncertainties in predicting the num- 
ber counts of galaxies. Sources of such uncertainties considered 
here include the time scale of star formation (§4.2. 1), the galaxy 
size (§4.2.2), and the adopted SN feedback (§4.2.3). 

4.2.1. Star Formation Time-Scale 

Here we evaluate the effects of changing the time-scale of 
star formation from CSF (equation 1) to DSF (equation 2). Fig. 
A7 shows predicted number counts for the corresponding A- 
dominated LC and LD models. The difference between these 
models is apparently small but its magnitude-dependence is dif- 
ferent among the predictions in the UBVI bands. In the case of 
longer wavelength such as the hu band, the number of faint 
galaxies in the LD model becomes larger than that in the LC 
model, because in the case of DSF more stars are formed at 
high redshift according to much a shorter (DSF) as compared 
with CSF. [In the particular case of the EdS universe, not shown 
here, the DSF gives r*(DSF) oc (1+z)" 3 / 2 , so that the present 
t* (DSF) is more than 10 times longer than that at z <~ 5.] How- 
ever, in the case of shorter wavelengths, particularly the t/300 
band, because the apparent luminosity from galaxies is domi- 
nated by instantaneous SF, and because many stars have already 
been formed at higher redshift, the number of faint galaxies 
in the LD model becomes slightly smaller than that in the LC 
model. 

The difference between the LC and LD models is most 
prominent when we consider the redshift distribution of faint 
galaxies. Fig. A8 shows such predictions for the LC and LD 
models. Clearly the LD model predicts too many high-z galax- 
ies to be consistent with observation. We may be able to remedy 
this defect by imposing a more efficient internal dust absorption 
in order to decrease the number of high-z galaxies. However, 
this manipulation simply decreases the number of faint galax- 
ies below that observed in the f/3oo- Thus, we suggest that the 
z-dependence of SF should be negligible and the CSF is a rea- 
sonable option in the framework of our SAM. 

4.2.2. Galaxy Size 

In this paper, it is assumed according to the usual SAM anal- 
ysis that the physical mechanism for determining the galaxy 
size is the conservation of specific angular momentum for spi- 
ral galaxies and the virial theorem for the baryon component of 
elliptical galaxies. While this assumption is considered to be 
reasonable, it is important to note the uncertainties in the nor- 
malization of the above relations. In this section we consider 
two LC variants and compare them with the original LC model. 
Figs. A9 and A10 present theoretical predictions from these 
three LC models. We rather arbitrarily decrease and increase 
the original galaxy size by a factor of 2, and refer to these vari- 
ants as "high surface brightness" and "low surface brightness," 
respectively. Note that the la scatter in the observational data 
is about a factor of 1.7, so the range of changing the radius 
is sufficient large to check the uncertainty in the galaxy size. 
Given the threshold of surface brightness for detection used in 
the HST observation, model galaxies with "high surface bright- 
ness" are more easily detected, while those with "low surface 
brightness" remain undetected. The number counts of these two 
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variants differ only by a factor of 1.5 at B450 <~ 28, as well as 
at the corresponding magnitudes in the other bands. We found 
that the uncertainties from changing the galaxy size in the SC 
and OC models are almost the same as those in the LC model. 

4.2.3. Supernova Feedback 

The SN feedback-related parameters essentially determine 
the resulting shape of the local luminosity function of galaxies 
(§3.1). Therefore, changing these parameters corresponds to 
changing the number density of galaxies in the local universe. 
In this section, we constrain these parameters against the ESP 
luminosity function (Zucca et al. 1997) which gives the steep- 
est faint-end slope among other observations. This is done by 
weakening the strength of SN feedback in three CSF variants 
such as SCw, OCw, and LCw. Other parameters such as the SF 
rate-related are determined by the same way in §3. All the pa- 
rameters for these models are tabulated in Table A2, and their 
B-band luminosity functions are shown in Fig. All. 

Straightforward calculations predict the number counts of 
galaxies which exceed the bright counts at B <~ 20, and this ex- 
cess alone may invalidate the extreme assumption for faint-end 
slope of the local luminosity function. However, considering 
the possible existence of systematic uncertainties that could af- 
fect current observations of the local luminosity function, we 
adjust the normalization in such a way as to reproduce the ob- 
served bright counts. Then we examine the effect of adopting 
the steepest ESP slope instead of our standard choice. Fig. A12 
shows the number counts of galaxies for the SCw, OCw, and 
LCw models. It is evident from this figure that even the steep- 
est ESP slope does not save the standard CDM universe, which 
confirms the claim by TY00. 

5. SUMMARY AND DISCUSSION 

We have calculated the number counts of faint galaxies in the 
framework of a SAM for three cosmological models of the stan- 
dard CDM (EdS) universe, a low-density open universe, and a 
low-density flat universe with nonzero A. The novelty of our 
SAM analysis is that theoretical predictions are made by fully 
taking into account the selection effects from the cosmological 
dimming of surface brightness of galaxies and also from the 
absorption of visible light by internal dust and intergalactic H I 
clouds. 

Comparison of theoretical predictions with the observed 
number counts and photometric redshift distribution of HST 
galaxies, as well as other ground-based observations, indicates 
that the standard CDM is ruled out and a A-dominated flat uni- 
verse is most favorable, while a low-density open universe is 



marginally favored. This result is in sharp contrast with pre- 
vious SAM analyses on galaxy number counts where many 
of the conceivable selection effects in faint observations have 
been ignored. It is only recently that the SAM analyses have 
included the effects of internal dust absorption (Somerville & 
Primack 1999; Cole et al. 2000a), and intergalactic H I absorp- 
tion (Baugh et al. 1998), and the isophotal selection effect in a 
very simple way (Cole et al. 2000a). However, as stressed by 
TY00, any predictions based on number count analyses will be 
seriously compromised unless all of the selection effects con- 
sidered in this paper are taken into account simultaneously. 

Based on a hierarchical clustering scenario in the CDM 
universe, the SAM naturally involves the merger-driven num- 
ber evolution of galaxies, which has been introduced only 
phenomenologically in traditional models of galaxy evolution 
(Yoshii 1993; Yoshii & Peterson 1995; TY00). The fact that we 
have essentially reached the same conclusion from phenomeno- 
logical approach confirms that previous simple prescriptions of 
the number evolution of galaxies are still useful in studying the 
global evolution of faint galaxies. It should be noted that in 
TY00 the adopted number evolution law which reproduces the 
observational data is oc (1 +z) andL* oc (l+z)~\ and is con- 
sistent with the observational estimate of the merger rate of the 
Canada-France redshift survey (CFRS) galaxies (Le Fevre et al. 
2000). 

The basic ingredients in the SAM analysis include the SF 
process and SN feedback. Although the constant and dynami- 
cal SFs can equally reproduce the local luminosity function by 
adjustment of their free parameters, we find that the dynamical 
SF predicts the formation of too many high-z galaxies to be con- 
sistent with the photometric redshift distribution of faint HST 
galaxies. Thus, our SAM analysis prefers the constant SF, as is 
also supported from other recent SAM analyses on the forma- 
tion of quasars (Kauffmann & Haehnelt 2000) and the evolution 
of damped Ly-a systems (Somerville, Primack & Faber 2001). 

The SN feedback, associated with the virial equilibrium for 
the baryonic component, controls the resulting size versus lumi- 
nosity relation of elliptical galaxies. However, our SAM analy- 
sis is unable to reproduce the observed relation bifurcating into 
the giant-dwarf and giant-compact sequences; our brighter el- 
lipticals reside on the bright portion of giant-dwarf sequence 
and our fainter ellipticals on the faint portion of giant-compact 
sequence. Although this limitation of our SAM analysis is 
found to hardly affect the conclusion in this paper, it is of urgent 
importance for the SAM to be equipped with some mechanism 
enabling the bifurcation of early-type galaxies into two distinct 
sequences as observed. 



APPENDIX 

A. BULGE SIZE AND SUPERNOVA FEEDBACK 

In §4.2.2 we mentioned that the intrinsic size of early-type galaxies is related to the strength of SN feedback. In §2. 1 describing 
SN re-heating the baryonic mass in a halo is given by Mb < SlbM H / (1 +/?), where (3 = t^/t^^ and M H is the halo mass. The left- 
and right-hand sides are not equal to each other, because this relation depends on the merging history of halos. The spherical collapse 
model gives Mh oc V^, ignoring the dependence on the formation redshift. Thus, we obtain r e oc Mb/Vj rc oc V C irc/(l +/?)• Assuming 
a constant ratio of baryonic mass Mb relative to luminosity L, two limiting cases of [3 <C 1 and /? 3> 1 give 

--^—M + const. for (3 <C 1 

7 fl+«hot) onst for /3>>1 < A1 > 



log r e = < 



2.5(3 + a hot )' 

where M is the absolute magnitude. Fig. A13 shows the size of early-type galaxies for the SC model (upper thick curve) and the 
LC model (lower thick curve). The vertical scale is chosen arbitrarily to avoid overlapping. Solid lines are for (3 <C 1 in equation 
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(Al). Dot-dashed and dashed lines are for (3 3> 1 with ahot = 5.5 and 2.5, respectively (see Table Al). Although we fixed the 
baryon fraction, the baryonic mass-to-light ratio, and the formation redshift of halos, theoretical results (thick curves) can readily 
fit the scaling relations. This indicates that the SN feedback process essentially determines the size of early-type galaxies. In other 
words, we can constrain the SN feedback from reproducing their observed size versus luminosity relation. It should be noted that the 
importance of SN feedback in galaxy evolution has also been highlighted in other recent SAM analyses of the color versus magnitude 
relation of early-type galaxies (Kauffmann & Chariot 1998; Nagashima & Gouda 1999). 

We would like to thank T. C. Beers for his critical reading of the manuscript. This work has been supported in part by the Grant- 
in-Aid for the Center-of-Excellence research (07CE2002) and for the Scientific Research Funds (10640229 and 12047233) of the 
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FIG. A 1 . — Extinction in units of magnitude by intervening intergalactic H I 
clouds for four HDF passband filters. The model optical depth is taken from 
Yoshii & Peterson (1994). 
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FIG. A2. — Local luminosity functions in the (a) B band and (b) K band. The 
solid, dot-dashed, dashed, and dotted lines indicate the SC, OC, LC and LD 
models, respectively. Symbols with errorbars in (a) indicate the observational 
data from APM (Loveday et al. 1992, filled squares), ESP (Zucca et al. 1997, 
open triangles), Durham/UKST (Ratcliffe et al. 1998, open diamonds), 2dF 
(Folkes et al. 1999, open circles), and SDSS (Blanton et al. 2000, crosses). 
Note that the SDSS luminosity function shown here is that with the same de- 
tection limit as employed in the 2dF survey. Symbols in (b) indicate the data 
from Mobasher et al. (1993, open squares), Gardner et al. (1997, open stars), 
and 2MASS (Cole et al. 2000b, filled triangles). 
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FIG. A3. — Cold gas mass relative to B-band luminosity of spiral galax- 
ies. Solid curves in the four panels from top to bottom show the SC, OC, LC 
and LD models, respectively. Errorbars denote lcr scatters. Filled diamonds 
indicate the observational data for atomic neutral hydrogen taken from Hucht- 
meier & Richter (1988). In the models, the cold gas consists of all species 
of elements, therefore its mass is multiplied by 0.75, i.e., Mm = 0.75M co id, 
which corresponds to the fraction of hydrogen. Because the observational data 
denote only atomic hydrogen, they should be interpreted as lower limits of the 
ratio. 
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FIG. A4. — Galaxy sizes. The thick solid lines and the thick dashed lines show theoretical results for effective radii of late-type and early-type galaxies, 
respectively. Errorbars denote \a scatters. The thin dotted lines, the thin dashed lines, and the thin dot-dashed lines show the observed best-fit relations given by 
TY00 for spirals, dwarf ellipticals and compact ellipticals, respectively. The original data are given by Impey et al. (1996) for spirals and Bender et al. (1992) for 
ellipticals. 
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FIG. A5. — Number-magnitude relations for various cosmological models. The solid, dot-dashed, and dashed lines indicate SC, OC, and LC, respectively. The 
thick lines denote models including the selection effects and the absorption effect by intergalactic H I clouds, while the thin lines denote those without the effects. 
The symbols indicate observational data. 
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redshift distribution of the HDF. 
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FIG. A9. — Number-magnitude relations. The solid lines indicate the original LC model. The thick and thin lines denote models with and without the selection 
effects. The dashed lines and the dot-dashed lines indicate the high and low surface brightness models, respectively. 
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FIG. A10. — Redshift distributions. The thick solid lines indicate the LC model. The dashed and dot-dashed lines denote the high and low surface brightness 
models for the model LC, respectively. The histograms indicate the photometric redshift distribution of the HDF. 
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FIG . All. — Local luminosity function in the B-band. The solid, dot-dashed and dashed lines indicate the SCw, OCw and LCw models, respectively. The types 
of symbols denote observations and are the same as Fig. A2a. 
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FIG. A13. — Galaxy sizes in arbitrary unit. The upper and lower thick lines indicate effective radius of early-type galaxies in the SC and LC models, respectively. 
The thin solid lines indicate the case of <C 1 and the thin dot-dashed line and the thin dashed line indicate the case of /3 2> 1 with Qhoi = 5.5 and 2.5, respectively. 
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Model Parameters 





cosmological parameters 




astrophysical parameters 








CDM Model 
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n A 
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°"8 


Vhot (km s ') 




< (Gyr) 




/bulge 


fb 


T 


SC 


i 





0.5 


0.67 


320 


5.5 
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-3.5 


0.2 
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1. 


OC 


0.3 





0.6 
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220 
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-3 


0.5 


0.5 


1.5 


LC 


0.3 


0.7 


0.7 


1 


280 


2.5 


1.5 


-2 


0.5 


0.5 


1.5 


LD 


0.3 


0.7 


0.7 
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280 


2.5 
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-2 


0.5 


0.5 


1.5 
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Parameters for Weak Feedback Models 












cosmological parameters 




as 


trophysical parameters 








CDM Model 
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n A 
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°"8 


V hot (kms- 1 ) 


"hot 


< (Gyr) 


a* 
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fb 
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SCw 


i 





0.5 


0.67 


300 


4 
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-2.5 
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OCw 


0.3 





0.6 
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160 


2.5 


1.5 


-1.5 


0.5 


0.5 


1.8 


LCw 


0.3 


0.7 


0.7 
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160 
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1.5 


-1.5 


0.5 


0.5 


1.5 



